Perturbation of magnetostatic modes observed by FMRFM 



in 
o 
o 

(N 
5h , 

< 

(N 
(N 



R. Urban 1 , A. Putilin 1 , P.E. Wigen 1 - 2 , S.-H. Liou 3 , M.C. Cross 1 , P.C. Hammel 2 , and MX. Roukes 1 
1 Condensed Matter Physics, California Institute of Technology, MC 114-36, Pasadena, CA 91125, USA 
2 Department of Physics, Ohio State University, Columbus, OH 43210, USA 
3 Department of Physics and Astronomy, University of Nebraska, Lincoln, NE 68588, USA 

(Dated: February 2, 2008) 

Magnetostatic modes of Yttrium Iron Garnet (YIG) films are investigated by ferromagnetic reso- 
nance force microscopy (FMRFM). A thin film "probe" magnet at the tip of a compliant cantilever 
introduces a local inhomogeneity in the internal field of the YIG sample. This influences the shape of 
the sample's magnetostatic modes, thereby measurably perturbing the strength of the force coupled 
to the cantilever. We present a theoretical model that explains these observations; it shows that 
tip-induced variation of the internal field creates either a local "potential barrier" or "potential well" 
for the magnetostatic waves. The data and model together indicate that local magnetic imaging of 
ferromagnets is possible, even in the presence of long-range spin coupling, through the induction of 
localized magnetostatic modes predicted to arise from sufficiently strong tip fields. 
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In the past decade, the possibility of using magnetic 
degrees of freedom in electronic devices has attracted the 
attention of the semiconductor industry. There are nu- 
merous proposals to incorporate spin degrees of freedom 
into electronics creating spin-electronic or spintronic de- 
vices Q. These devices require a thorough understand- 
ing of the magnetization dynamics on the submicron 
scale. Magnetic resonance force microscopy (MRFM) 
has proven itself to be a three dimensional nondestructive 
imaging technique that can be applied to electron 0,0, Q 
and nuclear 0] resonance experiments. In MRFM a mi- 
cromagnet (tip magnet) is placed on a compliant can- 
tilever. The time variation of the amplitude of the pro- 
cessing spins will produce a force on the cantilever via the 
gradient field of the tip magnet. In addition the tip mag- 
net produces a perturbation in the magnetic field in the 
sample producing a small bowl shaped region in which 
the spins are at resonance, the resonance slice. This res- 
onance slice can be highly localized (on the order of the 
size of the tip magnet) thus allowing subsurface detection 
of spins with submicron resolution. A recent experiment 
reported the detection of a single spin in resonance 0. 

Ferromagnetic Resonance Force Microscopy (FM- 
RFM) H is a variation of MRFM that enables the 
three dimensional, high sensitivity and high resolution 
qualities of MRFM to be applied in the characterization 
of the dynamic magnetic properties of ferromagnetic thin 
films and multilayers at the micron scale. The tip field 
penetrates magnetic and nonmagnetic materials allow- 
ing one to investigate buried structures and interfaces. 
Recent publications have shown the ability of FMRFM 
to determine the nature and magnitude of terms con- 
tributing to the internal field 0, dispersion relations of 
magnetostatic modes E l B and relaxation processes 
in micron size samples 10] . Ferromagnetically coupled 
systems pose unique challenges for magnetic resonance 
imaging due to the strong coupling between the spins. 
The resulting magnetostatic/exchange resonance modes 
involve spins occupying the entire sample not just those 



within the resonance slice jfj. Therefore, local imaging 
capability of MRFM is lost in the case of ferromagnetic 
samples. 

In this Letter we report the first direct observation of 
the variation of the magnetostatic mode amplitudes due 
to the local perturbation of the internal field. The local- 
ized field is produced by the sharp magnetic tip when it is 
scanned across the surface of Yttrium Iron Garnet (YIG) 
samples. The experimental results are in excellent agree- 
ment with theoretical predictions. In addition, based on 
our theoretical model, we propose that it is possible to 
generate a highly localized mode and recover, in part, 
the advantages of local imaging capabilities of MRFM. 
These localized magnetic modes will provide an exciting 
opportunity to investigate static and dynamic magnetic 
properties of micron-size defects (imperfections) embed- 
ded in ferromagnetic samples without the interference of 
edge effects. 

The experiments were performed using an ambient 
FMRFM system in the perpendicular geometry 0. The 
sample is attached to a microstrip resonator having a 
resonance frequency of 7.4 GHz. The silicon cantilever's 
resonance frequency, the spring constants, and the ambi- 
ent Q-factor are 17.8 kHz, 0.4 N/m, and 90, respectively. 
The cantilever displacement is monitored by 850 nm fiber 
optic interferometer. The cantilever tip was coated with 
a 250 nm CoPt film (M s = 800 emu/cm 3 ) and annealed 
in an external magnetic field of 80 kOe oriented along the 
tip axis of the cantilever. This resulted in a coercive field 
of 10 kOe [llj- The magnetic tip generates a localized tip 
field, htip, which can be aligned either parallel (hf lp > 0) 
or antiparallel (h t i p < 0) to the external dc field, H ext - 
Based on a model of the magnetic tip, the strength of 
the tip field and the tip field gradient at the tip-sample 
separation z = 1 (an were estimated to be 80 G and 50 
G/^m for h t ip > 0. 

In the following section we present experimental data 
that confirm the importance of local probe-induced phe- 
nomena. Typical measured spectra obtained on the 



25 10 5 G 



h = 80 25 10 5 G 



1.0 

0.8 

0.6 

0.4 

0.2 

0.0 
1.0 

0.8 

0.6 

0.4 

0.2 

0.0 



- (a) z = 2 \im 

z = 10 um 

CO *- 

A'.jv\i.3V. A_. A (V . 






(b) - 

C 


I 1 I 1 I 1 

z = 3 urn - 

ST 





-300 -200 -100 

Sweep field [G] 

FIG. 1: Typical FMRFM spectra measured in the center 
(y = fim) of the 20 x 80 /im 2 sample for (a) ht\p > and 
(b) /itjp < 0, respectively. The dashed and the solid lines 
correspond to the tip-sample separation of 10 /im and 2 /im, 
respectively. Zero sweep field corresponds to the resonance 
field of the homogeneous mode w/7 + 4irM 3 — 4.4 kG. 



20 x 80 /im 2 YIG sample are shown in Fig. ^ The 
magnetostatic modes are labelled by (n x ,n y ), where n x 
and n y represent the number of half wavelengths along 
the sample width and length, respectively. In the pres- 
ence of a uniform bias field and a homogeneous driving 
rf field only modes with odd values of n x and n y have 
a net dipole moment that will couple to the rf field [l^ . 
However, if the symmetry of the magnetic field is broken 
by the presence of the localized tip field, the magneto- 
static modes having n x and/or n y even can have a net 
dipole moment and be excited. We now discuss the cases 
of parallel and antiparallel tip orientations. 

htip > 0: To investigate the effects of this perturba- 
tion quantitatively, we first discuss the case when the 
magnetic field of the tip magnet adds to the bias mag- 
netic field. For a weak tip field (h tlp ~ 0) applied at 
the center of the sample the mode amplitudes are ob- 
served to decrease monotonically with increasing n x and 
n y . Such a spectrum is shown by the dashed line in Fig. 

for z = 10 /im. At that distance the magnetic field 
generated by the tip magnet is on the order of 2 G. As 
the tip approaches the sample surface, the magnetic tip 
field increases. This results in a dramatic change of the 
measured spectra as illustrated by the solid line in Fig. 
[TJi. There are two main features to be noted: (i) the 
spectrum is shifted to a lower resonance field by (8 ± 1) 
G and (ii) the amplitude of modes (1,5) and (1,9) are 
enhanced relative to modes (1,3), (1,7), and (1,11). The 
resonance field shift is due to the spatially-averaged tip 
field which adds to the external field. Notice, that the 
principal mode (n x = n y = 1) is also split. We attribute 
the mode on the high field side to be a localized surface 
mode |l3j. 
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FIG. 2: (a) The mode intensity of the (1,3) and (1,5) mode as 
a function of the tip-sample separation. The tip was placed 
above the middle of the 20 x 80 fim 2 sample. The solid lines 
are a spline fit to data to guide the eye. (b) The theoretical 
prediction using a cylindrical magnetic tip. 



To demonstrate the subtle details of the effect of the in- 
creasing tip field, the force intensity for modes (1,3) and 
(1,5) is plotted in Fig. 0as a function of the tip-sample 
separation. As the tip field gradient increases with de- 
creasing distance the FMRFM signal increases. However, 
the intensity of the (1,3) mode decreases dramatically for 
z < 3 /im and at z = 2.5 /im its intensity becomes even 
lower than the intensity of the (1,5) mode. For z < 2 /im 
the intensity of the (1,5) mode decreases as well. These 
effects can be explained by the model proposed below, 
see Fig. EJd. 

htip < 0: By reversing the orientation of the applied 
magnetic field without reversing hf lp , the magnetic field 
of the tip subtracts from the homogeneous bias field. For 
a tip-sample separation below 2 /im, the modes (1,3), 
(1,7), and (1,11) are now stronger compared to modes 
(1,5) and (1,9), see Fig. ^p. This is exactly opposite to 
what is observed for hf lp > 0. Table [fl summarizes the in- 
tensity ratios of the subsequent resonant modes when the 
tip field is weak and the magnetostatic modes are practi- 
cally unperturbed, ha p « line and for strong tip fields, 
htip > and h t i p < 0. There are two notable features: 
(i) the values in the first row are lower/higher/lower than 
values for h tlp ~ while the values in the third row are 
higher/lower/higher than the values for hti p « 0. This 
reinforces the conclusions made earlier that for /i t ; p > 0, 
the modes (1,5) and (1,9) are enhanced, while for h tlp < 
the modes (1,3) and (1,7) are enhanced, (ii) The first and 
third columns are ascending while the second column is 
descending. 
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TABLE I: The intensity ratios for modes (1,3), (1,5), (1,7), 
and (1,9) in the presence of a strong positive tip field (h t i p > 
0), a weak tip field (hu p ~ 0), and a negative tip field (/iti p < 
0). Corresponding spectra are plotted in Fig. Q 
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FIG. 3: The mode intensity for the (1,2), (1,3), and (1,5) 
modes as the tip is scanned along the long axis of the 20 x 
80 /im 2 sample. The external dc field and the probe field are 
parallel. The tip-sample separation z w 3 /im. (a) Measured 
response, (b) Calculated amplitude. 



The influence of the perturbation field of the probe 
magnet is further demonstrated in lateral scans taken 
along the long axis of the sample at z s=s 3 /im. Fig. |3k 
shows the experimental data of the force amplitudes for 
the (l,n y ) modes as a function of position. Three inter- 
esting features are observed: (i) As the tip is moved from 
the center of the sample, the field of the tip magnet has 
broken the even symmetry of the internal field. As a re- 
sult the "hidden" (1,2) mode is excited, (ii) For the (1,3) 
mode, the intensity at y — /im is strongly suppressed 
while at the position of the next maxima in the magne- 
tostatic mode at y = ±26 /jm the intensity is enhanced, 
(iii) Mode (1,5) shows very little variation as a function 
of position of the tip. Therefore, the lateral resolution 
can be estimated to be approximately 20 /im. This is 
in good agreement with the theoretical predictions (Fig. 

Eb). 

These unusual characteristics are understood within a 
theoretical model which is based on a linearized Landau- 
Lifshitz- Gilbert (LLG) equation of motion that include 
a spatially dependent tip field and Maxwell's equations. 
The external dc field is parallel to the sample normal and 
the driving rf field, h r { , lies in the plane of the sample and 
is considered to be homogeneous over the entire sample. 
The cantilever tip field, hf lp (r), is considered to be either 
parallel, h tlv > 0, or antiparallel, h tlv < 0, to the external 
field H cxt . 

The linearized LLG equation of motion and Maxwell's 
equations neglecting displacement currents can be writ- 
ten as (in cgs units) 



dM tt . a 

-^ 7 [MxH cffl + - 
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V x H = 
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V • (H + 4ttM) = 
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where 7 (= 1.84 x 10 7 Oe _1 s _1 ) is the absolute value of 
gyromagnetic ratio, M s (= 140 erg/cm 3 ) is the saturation 
magnetization, and a (= 0.0045) is the dimensionless 
Gilbert damping parameter. The internal effective field, 
if off , is given by 

H cfi = (h x + h lU h y , flext - 4ttM s + h tip (r)) . (3) 

Eq. (JTJ is solved in a small angle approximation; \m\ <C 
M s , where m represents the transverse rf component of 
the magnetization vector. The boundary conditions are 
assumed to be 



= m „ = 



y=±L y /2 



= 0. 



(4) 



The unique feature in the model is treating the tip field 
as a non-local perturbation and determining its effect on 
the magnetostatic modes. Using Eq. |2 the rf magnetic 
field h is expressed in terms of the magnetization m(r) 
as a linear functional h[m(r)] resulting in the non-local 
relationship 



Km(r) + h (r)m(r) = Afim(r) + h T { . 

.7 
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ho(r) 



(5) 
(6) 

(7) 



The operator K represents a linear integral transforma- 
tion defined by 



^4( r ) = w fc4( r ) 



(8) 



where represents the unperturbed magnetostatic 
modes given by 

f k (x,y) = sin {k x (x + L x / 2)} sin {k y (y + L y /2)} (9) 

and ujfe is their dispersion relation given by the Damon- 
Eshbach theory [14|. The wave- vectors k are chosen to 
satisfy the "pinned" boundary conditions, Eq. ^ The 
new basis of the magnetostatic modes m n (r) can be 
found by solving the homogeneous equations 



K + h (r) m n (r) = Att n m n (r) . 



(10) 



The FMRFM signal is proportional to the force acting 
on the cantilever which can be written as 



FMRFM 



M z (r)-y z h tip (r)dxdy, (11) 



where M z (r) — M s — \m\ 2 /2M s is the longitudinal com- 
ponent of magnetization, d is the thickness of the film 
and A represents the area of the sample. A detailed dis- 
cussion of the theory will be presented elsewhere . 

Eq. 1101 offers an interesting interpretation; It is sim- 
ilar to the time-independent Schrodinger equation with 
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FIG. 4: The calculated magnetostatic modes of the 20 x 80 
jum 2 sample with the tip located at the center of the film. 
The unperturbed (dashed) and perturbed (solid) modes are 
plotted for a potential barrier htip = +60 G. Notice, that the 
wavelength of the mode is increased in the perturbed region 
(shaded), (a) For the (1,5) mode this region increases the 
net dipole moment and therefore increases the FMRFM sig- 
nal, (b) For (1,3) mode this region decreases the net dipole 
moment and therefore reduces the FMRFM signal. 



the kinetic energy-like term given by lo^, in Eq. |SJ The 
normalized tip field h (r) in Eq. [7| plays the role of a 
potential energy. When the tip field is parallel to the 
external field, ho(r) > 0, it creates a potential barrier 
for the magnetostatic modes, while ho(r) < yields a 
potential well. This analogy provides qualitative insight 
into how the magnetostatic modes are modified by the 
tip magnetic field. 

fotip > 0: The effect of the perturbation field when 
the tip field is parallel to the bias field is to increase the 
wavelength of the magnetostatic mode in the region near 
the tip. Fig. 0] shows this effect when the tip is located 
at the center of the film. If m is positive (i.e., when m 
is parallel to the net dipole moment of the mode) in the 
region near the tip as in the case of mode (1,5) shown 
in Fig. the net dipole moment of the mode will be 
increased and therefore the FMRFM signal is enhanced. 
On the other hand, if m is negative (i.e., when m is 
antiparallel to the net dipole moment of the mode) in 
the region of the tip field as shown for mode (1,3) in Fig. 
|5>, the net dipole moment will be decreased and therefore 
the FMRFM signal for that mode is suppressed. 

As the tip approaches the surface of the sample the tip 
field gradient increases and therefore the FMRFM signal 
increases, see Fig. [21 However, when the field separa- 
tion of the magnetostatic mode 47rM s AQ (see Eq. [(jj) is 
smaller than the tip field /itip, the mode becomes evanes- 
cent. This decreases the amplitude of the magnetostatic 
wave in the vicinity of the tip magnet resulting in a re- 
duced FMRFM signal even though the gradient of the 
field is increased. This is demonstrated in Fig. 



When the tip magnet is scanned across the sample, the 
transverse component of the magnetization m changes 
sign at the next lobe position; At y = /jm, the trans- 
verse component of m for the (1,3) mode is negative 
and therefore the intensity of this mode is strongly sup- 
pressed. On the other hand, at y = ±26 /im, m is pos- 
itive and, as observed in Fig. 4 this mode is strongly 
enhanced. 

^tip < 0: In the case of the potential well the results 
are exactly the opposite to those for a potential barrier 
(h t ip > 0). If m is negative in the region of the tip mag- 
net as is the case for mode (1,3), the wavelength of the 
wave in the region of the perturbing tip field is decreased 
and the coupling of the mode to the rf field is increased. 
Consequently the FMRFM intensity of the mode is en- 
hanced. On the other hand, if m is positive in the region 
of the tip field as for mode (1,5), the net dipole moment 
will be decreased and therefore the FMRFM signal for 
that mode is suppressed. 

Extending the Schrodinger equation analogy to the 
case of a deep narrow potential well, a highly localized 
magnetic mode can be excited |l5j . The properties of 
this localized mode would be independent of the sample 
size and therefore of any edge imperfections that might 
be introduced in the sample preparation. Scanning the 
local mode about the sample would allow one to ex- 
tract variations in the local magnetic properties of the 
ferromagnetic samples. This makes FMRFM a unique 
experimental method allowing the investigation of local 
magnetic properties of the ferromagnetic samples inde- 
pendent on the sample dimensions, shape and/or defects 
at the sample edge. Incorporation of advanced nano- 
magnetic probe tips should enable next-generation FM- 
RFM systems capable of imaging local magnetic proper- 
ties with sub-micron lateral resolution. 

In conclusion, it has been demonstrated that the mag- 
netic field of the tip magnet introduces a local inho- 
mogeneity in the internal field of a 20 x 80 /im 2 YIG 
film. This influences the shape of the samples mag- 
netostatic modes, thereby measurably perturbing the 
strength of the force coupling to the cantilever. Using 
a hard magnetic coating on the tip, we are able to in- 
vestigate the magnetostatic modes when the tip field is 
parallel (h t i p > 0) and antiparallel (h t i p < 0) to the ex- 
ternal dc field. The condition h t i p > can be visualized 
as a potential barrier in the Schrodinger-like equation 
analogy, and modes (1,5) and (1,9) are enhanced while 
modes (1,3) and (1,7) were suppressed. In contrast with 
this is the case for /i t i p < corresponding to a poten- 
tial well in the Schrodinger-like equation analogy. Here 
the observed effect is exactly the opposite to the case 
for htip > 0. These results are in excellent agreement 
with the proposed theory indicating, that an inhomo- 
geneity in the internal field can also be successfully in- 
cluded into FMRFM experiments. Furthermore, highly 
localized magnetostatic modes are predicted which are 
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suitable for local imaging of ferromagnetic samples. 
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